Interplay between Quantum Size Effect and Strain Effect on Growth of Nanoscale 

Metal Thin Film 
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We develop a theoretical framework to investigate the interplay 
and strain effect on the stability of metal nanofilms. The QSE 
coupled through the concept of "quantum electronic stress. First- 
quantum oscillations in the surface stress of metal nanofilms as 
adds extrinsically additional strain-coupled quantum oscillations 
nanofilms. Our theory enables a quantitative estimation of the 
samples, and suggests strain be an important factor contributin; 
existing theories and experiments. 
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When the thickness of a metal film is reduced to the 
range of electron Fermi wavelength, quantum confine- 
ment becomes prominent to form discrete quantum well 
states, giving rise to various manifestations of quantum 
size effect (QSE) [lf{ . In particular, the QSE has been 
shown to be a dominant factor in the growth of metal 
nanofilms on semiconductor substrates [lH5j in the so- 
called electronic growth regime [2|. On the other hand, 
the strain effect is ubiquitous in heteroepitaxial growth 
of semiconductor and metal thin films [g, [7[ . A few re- 
cent studies (il-H^I have considered both effects on metal 
thin film growth. One thermodynamic theory [i( stud- 
ied both effects on film stability, and two kinetic models 
Tol 11 1 assumed growth parameters to be dependent of 



island height and radius due to the QSE and strain ef- 
fects. However, majority studies have focused on one 
effect while neglecting the other, and those few studies 
which considered both effect have been generally limited 
to treat them as two independent additive effects. This 
is mostly because fundamentally no theory is available 
to assess how the QSE may change the stress state of 
the film, and conversely how strain may alter the QSE. 
Therefore, it is very important to establish a theoretical 
framework that underlies the QSE on surface stress that 
in turn unerlies the interplay between the QSE and strain 
effect. 

The Pb(lll) film grown on Si(lll) substrate has 
been extensively studied as a model system for QSE |3j- 
& EH 0] ■ The almost perfect matching between the Pb 
Fermi wavelength and its interlayer spacing in the (111) 
direction gives rise to two striking QSE features in Pb 
film: the odd-even oscillations and beating patterns ex- 



hibited in many properties such as surface energy and 
stability. These two main features have been agreed upon 
by all theoretical and experimental studies (JIB, IH 14 1. 
However, there remain some outstanding discrepancies. 
Oscillation patterns may vary slig htly from one experi- 



mental sample to another 13H16T]. First-principles cal 
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culations [4j predicted that the odd-even oscillations in 
surface energy essentially die out at a thickness of ^20 
monolayers (MLs), while experiments, in contrast, have 
seen the large oscillations sustain even beyond 30 MLs 
(l3l EH ■ One origin of the discrepancies was attributed 
to Pb/Si interface that causes a phase shift in the oscil- 
lation patterns |3j, but the strain effect has been mostly 
overlooked so far. 

Because of the large lattice mismatch, the Pb (lll)film 
tends to grow on Si(lll) substrate by adopting a 10-to-9 
epi relation to minimize interfacial misfit strain [It], EH ■ 
Even so, Pb film can still experience up to ±3% strain de- 
pending on the film orientation relative to Si surface [13] ■ 
The measurement of interlayer spacing by X-ray diffrac- 
tion 13j suggested that the strain in Pb film be small 
based on bulk Poisson ratio, but the actual amount of 
in-plane strain remains uncertain, because the ultrathin 
film may not follow the bulk Poisson ratio, especially in 
the presence of QSE that modifies the interlayer spacing. 
Overall, the strain effect has not been studied adequately 
in relation with the QSE, because of the lack of theory 
underlying their relationship and because the direct mea- 
surement of strain in the film is very difficult. 

In this letter, we develop a general theory underly- 
ing the fundamental relationship between the QSE and 
strain effect in the formulation of surface energy through 
the concept of "quantum electronic stress" [19j, i.e. the 
additional surface stress oscillations induced by the QSE. 
It allows us to theoretically study the interplay of these 
two effects on the stability of metal nanofilms by treating 
both effects on the same footing. Using first-principles 
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calculations, we reveal large quantum oscillations in the 
surface stress of Pb(lll) films as a function of thick- 
ness, which adds extrinsically additional strain-mediated 
quantum oscillations to surface energies of the strained 
Pb films. Our theory enables a quantitative estimation 
of the amount of strain in different experimental samples 
from the measured stability patterns. 

We first briefly introduce the concept of quantum elec- 
tronic stress that gives rise to the quantum oscillations 
of surface stress. Following density functional theory 
(DFT), the total energy functional of a solid is written 
as 

E[n(r),{R m }} = E e [n{r)]+E ext [n(r) 1 {R m }}+E I [{R m }} 

(1) 

E e [n(r)\ is the electronic energy functional of charge 
density n(r), including kinetic and electron-electron in- 
teraction energy, E ext [n(r), {i? m }] is the ion-electron in- 
teraction energy, Ei[{R m }\ is the ion-ion interaction en- 
ergy and {R m } are atomic coordinates. Consider a vari- 
ation of electron density from the ground-state n° as 
n* = n°+Sn in the absence of strain (i.e., without any lat- 
tice deformation), a general expression for lattice stress 
induced by such pure electronic perturbation/excitation 
can be derived as fijl 
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(2) 



which is called quantum electronic stress, [i is electron 
chemical potential and dfi/dsij is electron deformation 
potential. In a nanofilm of thickness d, QSE induces vari- 
ation of charge density and deformation potential along 
the surface normal z-direction. Then a special form of 
quantum electronic "surface" stress can be expressed as 



qe 1 f d[i 
°« = d J d^ Sn{z)dz - 



(3) 



We have performed DFT calculations to directly reveal 
quantum surface stress oscillations in Pb(lll) nanofihns. 
Our calculations are done using VASP code [20| based on 
density functional theory in plane-wave formalism. For 
all the freestanding Pb films and Pb film on Si substrate 
from 1~11 MLs, ultrasoft pseudopotential [2lj and gen- 
eralized gradient approximation are used with the Pb 
5d orbitals included as valence states. For thicker Pb 
film (12MLs and thicker) on Si substrate, PBE poten- 
tial [22I ] and generalized gradient approximation without 
5d orbitals are used to save time. All calculations use a 
plane-wave cutoff of 240eV to obtain good convergence 
for stresses which typically converge slower than total 
energy. The Pb film is modeled by a supercell slab with 
the strain-free film set at the theoretical lattice constant 
of 5.04A. The Si substrate was modeled using 6 layers 



140 - 



120 



> 
CD 

100 

b 



80 




a i +(J QE 











if 




jfjww — — 




60 

20 40 60 80 100 120 140 160 180 

d (ML) 

8 i u 400 

-320 

-240 

-160 q 
o 

-80 ~£ 
O 

-0 < 
--80 ^ 
--160 
--240 
-320 

2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 

d(ML) 

FIG. 1: Surface energy (squares) and surface stress (dots) 
of freestanding and Si-supported unstrained Pb(lll) film ob- 
tained from DFT calcualtions. The insets show schematics of 
film. 



of Si with the bottom two layers fixed at bulk positions 
and the bottom layer passivated with H. The slabs are 
separated by a vacuum thickness of ^20A in z-direction, 
sampled by a 20x20x1 mesh in k-space. 

Figure 1(a) shows the calculated surface stress (<Jo), as 
a function of film thickness (d) up to 130 MLs, of the 
freestanding strain- free Pb(lll) film. It is well-known 
that surface energy displays an oscillatory dependence on 
d [4|. What's new is that surface stress a displays also 
a strong oscillatory dependence on d, indicating a strong 
QSE on surface stress, as suggested by Eq. (3) above. 
The surface stress exhibits a similar oscillation pattern 
with film thickness as surface energy characterized 
by an odd-even oscillation superimposed by a beating 
pattern with a period of ^9 layers, but the oscillations 
of surface stress are out of phase with those of surface 
energy and the beating patterns are phase shifted too 
(not shown). 

where E c is the classical surface energy (bond break- 
ing energy) of a macroscopic thick film independent of 
film thickness, and E® is the quantum surface energy 
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due to the QSE, which is function of film thickness d. 
By definition, surface stress (a rank-2 tensor) can be ex- 
pressed as 

In general, we may also express the surface stress as 



(<9, 



(4) 



where a M is the mechanical surface stress of a macro- 
scopic thick film which we are familiar with, and a® E 
is the new oscillating quantum surface stress as a func- 
tion of film thickness d. The thickness dependence of the 
quantum surface stress can be related to the thickness 
dependence of charge density and electron deformation 
potential induced by QSE as shown in Eq. (3). As the 
film thickness increases, however, <jQ E will eventually di- 
minish and a will converge to a M as indicated in Fig. 
1(a) although we could't calculate thicker film beyond 
the 130 ML to show full convergence. 

The introduction of the quantum surface stress pro- 
vides a direct link between the QSE and strain effect on 
the surface energy and hence stability of thin films in the 
quantum regime. In particular, under a given strain e, 
the surface energy will have the following thickness de- 
pendence within linear elasticity 



E(e) = E {d) + A[a M + a^ E (d)} ■ e 



(5) 



Where the first term is the surface energy of a un- 
strained film (denoted by subscript "0") which has a 
thickness (d) dependence (quantum oscillations) due to 
the QSE alone, while the second term is the strain in- 
duced surface energy which adds extrinsically additional 
strain-coupled quantum oscillations to surface energy 
through oscillating quantum surface stress. Equation 
(5) enables a quantitative assessment of the interplay be- 
tween the QSE and strain effect on the stability of metal 
nanofilms. 

In experiments, Pb films are grown on semiconductor 
substrates, such as Si and Ge. Hence, in order to compare 
with experiments, we must also include the substrate and 
interfacial effects. Figure 1(b) shows the calculated 70 
and (To as a function of d ranging from 1 to 31 MLs of the 
strainfree Pb(lll) film on a Si substrate (To do so, the Si 
substrate is strained to match the Pb lattice Q). There 
are some interesting similarities and differences between 
the freestanding and substrate-supported Pb (111) film. 
In both cases, 70 and ao show an odd-even oscillation 
modulated by a 9-layer beating pattern; do displays a 
larger oscillation magnitude than 70. On the other hand, 
for the substrate-supported film, 70 and 00 contain not 
only the surface terms but also the interface terms, which 
makes their oscillation amplitude almost twice as large as 
that of the freestanding film. In addition, the presence 
of Si substrate causes a phase shift in 70 and <tq by ~1 
ML; the first node of the beating pattern appears at the 
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FIG. 2: Comparison of surface energy between model pre- 
diction and direct DFT calculation for a Pb film under 1% 
strain, showing excellent agreement. 



7ML and 6ML for the freestanding and Si-supported Pb 
film, respectively. 

To verify our theoretical framework, we calculated the 
surface energies of the 1% strained film as a function 
of thickness in comparison with the model predictions, 
as shown in Fig 2. We see that the model predictions 
agree very well with the direct first-principles results, 
validating our theory. Thus, using the first-principles 
calculated the surface energies and surface stresses of the 
" unstrained" film, we can apply our model to predict the 
surface energy (7) of the strained film with or without 
substrate support. 

Figure 3(a)and(b) shows the predicted surface en- 
ergy of the freestanding and Si-supported Pb(lll) films 
strained from -3% to 3%, respectively. Strain modi- 
fies the surface energy in two important ways. First, 
strain enhances the QSE by increasing the odd-even os- 
cillation magnitude in 7. This enhancement extends the 
QSE induced surface energy oscillations to much thicker 
films (the oscillation persisting beyond 30 ML with ^3% 
strain). So, strain provides one possible reason for the 
experimentally observed stability oscillations existing in 
much thicker films (^30 ML) |13j| than the previous the- 
oretical predictions (~20 ML) Second, because the 
quantum oscillations in surface stress and surface energy 
are phase shifted, large enough strain will change the os- 
cillation pattern (both the odd-even and beating pattern) 
of surface energy. This means that strain will alter the 
relative film stability of different thicknesses. For exam- 
ple, for the strain-free freestanding film, the 14ML film 
is stable and the 15ML is unstable; however, under 3% 
strain, the 14ML becomes unstable and the 15ML be- 
comes stable, as shown in Fig 3(a). 

Experimentally, the observed stability patterns of 
Pb(lll) films grown on Si(lll) from different groups 
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FIG. 3: Model predicted surface energies of Pb(lll) films 
under strain from -3% to 3%. (a)freestanding film and (b)Si- 
supported film. 
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FIG. 4: Comparison of relative surface energies of Pb(lll) 
film on Si substrate between the experiment (dots) [13, L3 
and model prediction with the fitted amount of strain 
(square). 



are in generally good agreement but with some subtle 
differences around the nodal points of thicknesses in the 
beating pattern [13M16I ] . The reason for such discrepancy 
remains unresolved, although some general argument has 
been made by attributing the discrepancy to nonspher- 
ical Fermi surface [23[ and substrate effect 0, H| . Here, 
we argue that the discrepancy is partly caused by the 
different amount of strain in different experimental sam- 
ples. Below, we apply our model to extract the amount 
of strain in some experimental samples by matching the 
predicted stability patterns to the experiments. 

Without strain, the calculated stability pattern from 
either freestanding or Si supported film agrees poorly 
with the experiment by Czoschke 13 1 and Zhang [14|, 
as seen by comparing Fig. 1 with Fig. 4. In particu- 



lar, both experimental results show large odd-even os- 
cillations from 5 to 8 ML (Fig. 4), while the theory 
shows little oscillation in this region (Fig. 1) which is 
in the vicinity of a nodal point of the beating pattern. 
To resolve this discrepancy, we apply the above theo- 
retical framework [Eq.(5)] to predict the stability pat- 
tern of "strained" Pb films on the Si substrate, using 
the calculated surface/interface energies and stresses of 
the unstrained film on the Si substrate. In fitting the 
experimental data, we assume a nonuniform strain dis- 
tribution in the film that decreases linearly with the in- 
creasing film thickness [24, 25 1, and then treat the strain 
and its decay rate as fitting parameters. We obtained 
very good fitting results by using a linear strain profile 
of 1.76% - (d - 5) X 0.068% for Czoschke's sample [l3| 
and 1.80% - (d-5) X 0.061% for Zhang's sample [lj], re- 
spectively, as shown in Fig. 4. Most noticeably, our model 
correctly predicted the large odd-even oscillations in the 
range of 5-8 ML as seen in the experiments. This is be- 
cause there is a large oscillation in the surface stress in 
this range (see Fig.l), which induces additional oscilla- 
tions in surface energy when strain is applied. The fitted 
strain are only slightly different in the two samples by 
~ 0.1%, in accordance with the overall agreement be- 
tween the two experimental patterns. Surprisingly, this 
small difference is enough to account for the subtle differ- 
ences in the two experimental patterns in the thickness 
range of 12-14 ML, 21-23 ML and 30-31 ML, all in the 
vicinity of nodal points. Overall, the strain is small, less 
than 2% initially, and decays with the increasing film 
thickness to less than 1% beyond 10 ML and diminishes 
around 30 ML. The average strain in a 30 ML film is ~ 
0.9%, within the range of general estimation 17 1. 

Recently, Miller, et al. has shown a fundamental phase 
relationship between the oscillations of surface energy 
and of work function that their beating patterns are al- 
ways offset by 1/4 of a period [26[ . We have shown that 
the strain can not only change the odd-even oscillations 
but also shift the phase of beating patterns of surface 
energy [27]. Applying the Miller's phase relation to the 
Si-supported Pb(lll) film by assuming that the inter- 
face shifts the work function and surface energy phase 
together [28| , we can fit the phase of surface energy beat- 
ing pattern to match (by an offset of 1/4 of a period) the 
experimental phase of work function pattern, such as the 
one measured by Qi et al.||29|, using strain as a fitting pa- 
rameter. We obtained the best fit with an average 0.75% 
strain for this particular film, as shown in Fig. 5. 

In conclusion, we have developed a theoretical frame- 
work to investigate the interplay between QSE and strain 
effect on the thermodynamic stability of metal nanofilms, 
through the introduction of a new concept of "quantum 
stress" [l3|. In the present case, the quantum stress rep- 
resents the additional surface stress induced by QSE. 
Broadly, our theoretical framework can be extended to 
investigate the interplay between QSE and strain ef- 
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FIG. 5: Comparison of experimental work function pattern 
[29| with DFT calculated surface energy pattern without 
strain and with model predicted surface energy pattern with 
0.75% strain. Note the 1/4 of a period of phase shift between 
the experimental data and model prediction. 



feet on kinetic growth properties, such as surface diffu- 
sion and step-edge barrier, where quantum " diffusional" 
stress [30|| induced by the QSE can be derived from first- 
principles to play the role of quantum surface stress here. 
Thus, our theory will be applicable to both thermody- 
namic and kinetic properties of nanoscale thin films when 
QSE and strain effects are prominent. 
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